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O ■ ABSTRACT 

(N : 

^ I We use the redshift evolution of the Lya forest mean transmitted flux at z > 2.5 

I to infer the evolution of the intensity of the ionizing background, using theoretical 

' predictions for the density distribution of the intergalactic gas giving rise to the Lya 
forest. The proper background intensity declines gradually with redshift, decreasing by 

^ ■ a factor ~ 3 between 2 = 3 and z = 5. 

*^ . The gravitational structure formation theory of the Lya forest implies that, given 

, the observed mean transmitted flux of F = 0.1 in the highest redshift quasar (with a 

\ mean Lya forest redshift z = 5.2), only about 3% of the Lya spectrum should have a 

Q \ flux higher than 50% of the continuum, and less than 0.5% of pixels should have flux 

' above 75% of the continuum, assuming a uniform ionizing background intensity. We 

■ show that this is consistent with the spectrum shown in Fan et al. . 

2 ■ Subject headings: cosmology: theory — intergalactic medium — large-scale structure of 

c/3 . universe — quasars: absorption lines 

a ■ 



X 



1. INTRODUCTION 

The theory of the Lya forest based on structure formation in Cold Dark Matter models has 
been highly successful in accounting for the observable properties of the Lya forest, in particular 
the flux distribution and the power spectrum (see Croft et al. 1999, McDonald et al. 2000a, and 
references therein). As new quasars at ever higher redshift are discovered, the comparison of these 
observable quantities with the theory should be extended. Because the theory predicts the evolution 
of the spatial distribution of intergalactic gas, it is now possible to infer the evolution of the ionizing 
background intensity from the observed evolution of the mean flux decrement. Furthermore, as the 
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epoch of reionization is approached, fluctuations in the background intensity due to the discreteness 
of the sources should increase, affecting the Lya forest flux distribution. 

In this paper, we use the observations of three quasars at 2; > 5 (Songaila et al. 1999, hereafter 
SHCM; Stern et al. 2000; Fan et al. 2000) to infer the evolution in the hydrogen ionization rate, and 
to test if the flux distribution is consistent with theoretical predictions. In §2 we present predictions 
of the flux distribution at high redshift, using numerical simulations whose resolution and box size 
are shown to be sufficient to have produced convergent results. In §3 we compare the predictions 
to the observations to measure the evolution of the ionization rate and search for differences in the 
predicted and observed flux distribution. 



2. NUMERICAL RESULTS FOR THE FLUX DISTRIBUTION 

If we ignore the effects of peculiar velocity of the gas and of thermal broadening, then the 
optical depth to Lya scattering, r, at any given point of an observed spectrum is determined 
by the density of gas at the corresponding point in space, pg^s = A (QbPcrit), and by the gas 
temperature, which we assume to follow the relation T = Tq A^-^ (see Hui k Gnedin 1997): 

Here, tq is a constant, H{z) ~ Hq ^H'^ (1 + z)^^'^ is the Hubble constant at redshift z, Hq is 
the present Hubble constant, h = i?o/(100kms~^/Mpc), $7^ is the fraction of the critical density 
in pressureless matter, r_i2(-2;) is the photoionization rate of hydrogen in units of 10^^^ s^^, and 
(3 = 2 — 0.7 (7— 1) (where we use a{T) oc T~^''^ for the recombination coefficient). The dependence of 
the distribution and average of the transmitted flux, F = exp(— r) (which is the observed quantity), 
on the cosmological model can be understood from equation (1). Thus, as long as the distribution 
of A and the A — T relation are sufficiently well understood, the observed mean transmitted flux 
of the Lya forest can be used to determine the parameter combination {ni,h^)^/H{z)/r^i2{z) in 
equation (1). 

Knowledge of the distribution of A and T needs to come from the theory, with the use of 
hydrodynamic simulations. This theory can be tested by checking its predictions for the detailed 
form of the distribution of F and its power spectrum against the observations, and this comparison 
has shown that the Cold Dark Matter model with a cosmological constant that is supported by a 
large body of observational evidence in galaxy clustering and CMB fluctuations (Primack 2000, and 
references therein) is in good agreement with the observed Lya forest (Ranch et al. 1997, Croft et al. 
1999, McDonald et al. 2000a, and references therein). At the same time, the gas temperature over 
the relevant range of densities has been measured from the absorption line shapes (Ricotti, Gnedin, 
& Shull 2000; Schaye et al. 2000; McDonald et al. 2000b), although the systematic differences 
between these studies are still substantial. 
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We focus in this paper on the flux distribution at high redshift. This particular property of 

the Lya forest is affected only by the low-density intcrgalactic gas, given the large mean flux 
decrement (the reason is that high density gas causes absorption only in highly saturated regions 
of the spectrum) . We first show that this distribution is predicted without significant uncertainties 
due to the finite resolution and box size of the simulations. We use the Hydro-PM (HPM) code of 
Gnedin & Hui (1998), comparing it first to the fully hydrodynamic simulation in Miralda-Escude 
et al. (1996) to test that the approximate nature of the Hydro-PM method does not result in any 
significant differences in the predicted flux distribution at high redshift. When we create simulated 
spectra, we assume a A — T relation given by Tq = 2 x 10"^ K and 7—1 = 0, which is consistent 
with the observational determinations Ricotti et al. (1999) and McDonald et al. (2000b) at z ~ 4, 
although Schaye et al. (2000) find a lower temperature (for the calculation of the pressure term in 
the HPM simulations we assume T = before reionization at 2; = 7, and then interpolate between 
To = 25000, 7 - 1 = 0.0 at z = 7 and To = 19000, 7 - 1 = 0.2 at z = 3.9; however, we find that 
the results are almost completely insensitive to these assumptions, probably because the pressure 
does not have much influence in the voids). The cosniological model is ACDM, with 0,^ = 0.4, 
Oft/i^ = 0.015 (note, however, that the results we give for r_i2 are rescaled to ri^/i^ = 0.02), 
Cl\ = 0.6, h = 0.65, and n = 0.95, where is the fraction of the critical density in cosmological 
constant and n is the usual logarithmic slope of the primordial power spectrum. Simulated spectra 
are constructed in which the transmitted flux is binned into 2 A pixels, to match the observed 
spectrum in Fan et al. (2000) (we use this pixel size throughout the paper). In Figure 1, we show 
the flux distribution predicted for the full-hydro simulation as the heavy solid line. The light {long- 
dashed, dotted, solid, dashed) lines show HPM simulations with box size L =(17.8, 8.9, 4.4, 4.4), 
and a total number of particles N = (256^, 128^, 128^, 256^), respectively. All the results are 
at z = 4, but using the value of r_i2 that gives a mean flux decrement from the hydrodynamic 
simulation of F = 0.1, which matches the observations at z ~ 5. We keep r_i2 (rather than F) 
constant for all the simulations, so that we can see the difference in the predicted F. 

The result for the fully hydrodynamic simulation and for the HPM simulation that has almost 
the same box size and resolution (L = 8.9 and = 256^, not plotted) are essentially identi- 
cal. This justifies using the HPM simulations to examine the effects of changing the box size and 
resolution and to compute the flux distribution for comparison with observations at different red- 
shifts. The two simulations with the smallest box size, with resolutions L/N^/^ = 35 Kpc and 
17 h~^ Kpc (solid and dashed lines, respectively), show that simulations with inter-particle separa- 
tion of 35 Kpc have already achieved the required resolution; however, degrading the resolution 
to L/Afi/3 = 69/i~iKpc causes a change in F of 0.02, similar to the observational error bars of 
F. Comparison of the long-dashed and dotted (or the heavy and light solid) lines shows that the 
effect of increasing the box size beyond 5h~^ Mpc is small. We conclude that simulations with 
boxes L = 8.9 Mpc and number of particles N = 256^ are sufficient for our purpose to have 
achieved convergence. We find similar convergence in the HPM simulations at z = 5.2. We also 
find that different random realizations of the initial conditions of the simulation of this box size 
and resolution change F by ±0.01 (small compared to the observational error). 
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Fig. 1. — Cumulative distribution of transmitted flux. Heavy solid line: full hydrodynamic simula- 
tion. Thin lines: Hydro-PM simulations, with box sizes and particle numbers indicated. 
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3. RESULTS 

We now examine the flux distribution at high redshift predicted from the simulations to do 
a direct comparison to the observational constraints. We use the HPM simulation with L = 
8.9 Mpc and N = 256^, and create spectra binned into 2 A pixels. The curves in Figure 2(a) 
show the cumulative flux distribution at z = 5.2, for ten values of r_i2, corresponding to different 
values of F, indicated in the Figure. This can be compared to the spectrum in Figure 3 of Fan et 
al. (2000), from a quasar at z = 5.8 (where the mean redshift of the observed Lya forest, in the 
wavelength range 7140-7956 A used by Fan et al. 2000 to compute F, is 5.2). In order to match 
the observed F = 0.10 it 0.02, we need r_i2 ^ 0.17 ± 0.03, for the Lya forest model described in 
§2. 

Remarkably, there should be essentially no pixels in the spectrum at this high redshift where 
the flux comes close to the continuum level. Only 3% of the spectrum should have a flux higher than 
50% of the continuum, and less than 0.5% should have a flux greater than 75% of the continuum 
[see the inset of Fig. 2(a)]. In comparison. Fan et al. (2000) identify the minimum absorption in 
their observed spectrum to be r 0.4 (although they estimate that continuum uncertainty and 
noise make this number uncertain by about ±0.1). Absorption of r < 0.4 ov F < 0.67 occurs in 
approximately 1% of the pixels in our simulated spectra, so we predict about three such pixels 
to be found in the wavelength range 7300-7956 A examined by Fan et al. (2000). Therefore, 
the minimum absorption present in this spectrum is consistent with the predictions of the simple 
Lya forest model in our simulations. Wc find that decreasing erg by ~ 15% to match the results 
of McDonald et al. (2000a) increases the inferred value of r_i2 by ~ 30%, but does not affect 
significantly our conclusions about the flux distribution (once the observed F is matched). 

Next, the same lines are plotted in Figure 2(b) at redshift z = 4.5, which is about the mean 
redshift of the Lya forest in the z = 5.0 quasar spectrum that is shown in Figure 2 of SHCM. The 
mean transmitted flux in this spectrum is F = 0.25, which requires r_i2 ~ 0.21 (no error bar was 
given in SHCM for the value of F). According to SHCM, an upper bound on the minimum optical 
depth in their spectrum is r < 0.1. They identify this region as being produced by a void with 
density A ~ 0.2, to derive a lower limit to the background intensity of r_i2 > 2.8, for Sl^/i^ = 0.019. 

Our value of r_i2, inferred by matching F, is much smaller for the following reasons: Factors 
of ~ (20000/5000)-°-'^, (0.4/1.0)1/2^ and 0.65/0.5 are obtained from our different choices of Tq, Qrn, 
and h, respectively (note: we assume SHCM used h = 0.5, but their paper is not clear on this point). 
Our choices place more conservative lower limits on r_i2, and are more consistent with independent 
observations. The remaining factor 2.4 must be traced to the identification by SHCM of observed 
r = 0.1 with underdensity A = 0.2. In our simulations, A < 0.2 occurs in 1.7% of pixels in real 
space, and this probability is increased when the expansion of voids in redshift space is considered. 
Consequently, the minimum absorption in the 200A region considered by SHCM probably arises 
from a void with A < 0.2. Furthermore, the peculiar expansion of voids reduces the optical depth at 
fixed density. After the inclusion of all these effects through our simulations, and using r_i2 = 0.21, 
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Fig. 2. — Cumulative flux distributions for (a) z = 5.2 and (b) z = 4.5. Listed values of r_i2 and 
F apply to the lines from top to bottom. The inset box in (a) is an expanded version of the upper 
right hand corner of the full figure. 
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we predict ~ 0.4 regions with r < 0.1 in the wavelength range 6800 - 7000 A identified by SHCM. 

Lastly, the minimum absorption statistic is biased toward regions with strong noise contribution, 
so the value r < 0.1 is probably not conservative (our simulations give a completely acceptable 1% 
probability for r < 0.15). 

Finally, the quasar in Stern et al. (2000) at z = 5.5 has mean transmitted flux F = 0.10 ±0.02. 
The ionization rate required to match this F is r_i2 = 0.13 ± 0.03, at a mean redshift z = 4.93. 
Their spectrum as shown has only ~ 15 A resolution so we do not interpret it further. 

The evolution of the background intensity r_i2 derived from the observed mean flux decrement, 
for rtjjh^ = 0.02, is shown in Figure 3, using the values from the three quasars discussed here and 
the lower redshift results from McDonald et al. (2000a) (updated by McDonald et al. 2000b using 
the measured IGM temperatures). 

4. DISCUSSION 

There is a clear reduction of the proper intensity of the ionizing background with redshift, by 
a factor ~ 3 from z = 3 to z = 5. This reduction was noticed previously by Ranch et al. (1997) 
up to z = 4. With the new high redshift quasars, we find a continuation of this gradual decrease 
of r, but without a dramatic drop that would be the sign of approaching the reionization epoch. 
The consistency of the fraction of the Lya spectrum with low absorption with the observations also 
suggests that the background intensity fluctuations are not yet extreme at z = 5. 

The evolution of the background intensity is related to the evolution of the comoving emissivity, 
e, and the proper mean free path A, by e (1 + z)^ oc T \~^, if A is small compared to the cosmic 
horizon. Assuming that the rate of incidence of Lyman limit systems evolves as (1 + z)^-^ (Storrie- 
Lombardi et al. 1994), the mean free path of ionizing photons is A oc (1 + z)~'^, and therefore the 
comoving emissivity changes between z = 3 and z = 5 by a factor (1/3) (6/4)"^^^ = 0.5, indicating 
only a small reduction, which is clearly different from the reduction by a factor 0.14 predicted by 
the evolution of the comoving density of luminous quasars (Schmidt, Schneider, & Gunn 1995) and 
implies, therefore, that galaxies or other sources are dominating the emission of ionizing photons 
at 2; = 5. 

We thank Nick Gnedin for providing a copy of his Hydro-PM code and for helpful comments 
on the manuscript, Renyue Cen for the output from his hydrodynamic simulation, and Adam Steed 
for helpful comments on the manuscript. 

REFERENCES 

Croft, R. A. C, Weinberg, D. H., Pettini, M., Hernquist, L., & Katz, N. 1999, ApJ, 520, 1 
Fan, X. et al. 2000, AJ, submitted (astro-ph/0005414) 



-9- 




z 



Fig. 3. — Evolution of observed r_i2. Computed for Qsh^ = 0.02, and the ACDM model discussed 
in the text (changing Qb, for example, will shift all the points by the same factor). Uncertainties 
related to the quasar continuum and the gas temperature may introduce some uncertainty beyond 
the given error bars, especially for the 3 highest z points. 
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